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Abstract 

The radiative B c — * luj decay is investigated in the Standard Model in framework 
of the light cone QCD. The transition form factors and decay width are calculated. 
A comparison of light cone QCD and constituent quark model predictions on these 
quantities is presented. 
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1 Introduction 



The experimental and theoretical investigations of the heavy flavored hadrons constitutes 
one of the main research area in particle physics. This is due their outstanding role in the 
precise determination of the fundamental parameters of the Standard Model (SM), such 
as Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, leptonic decay constants etc., 
and for deeper understanding of the dynamics of QCD. In this sense, B c mesons occupy 
an exceptional place. Since B c mesons contain two heavy quarks, their decay channels are 
very rich compared to that of B q , (q = u, d, s) mesons. Moreover, QCD predictions on 
B c meson decays are more reliable since for the heavy quarks perturbation theory works 
quite well. Because of the above-mentioned reasons, the investigation of the properties of 
B c mesons receives special attention. Production and different decay channels of B c mesons 
are widely discussed in the current literature (for a review, see 0). In [fj, [§, the number of 
B c mesons that will be produced in LHC, is estimated to be ~ 2 x 10 s . This is clearly an 
indication of the real possibility of an experimental investigation of the properties of the B c 
mesons at LHC. In addition to that, the B c meson decay channels can bring about some 
background contribution to the B ± meson decays, with the same final states [|J, which is 
another reason that makes the precise study of the decay channels of B c mesons worthwhile. 

The pure leptonic decays of B c mesons are the simplest among all decays. In principle 
the pure leptonic decay B~ — > tv can be used for the determination of the leptonic decay 
constant fg c , which is one of the fundamental parameters of the hadron physics. But 
this type of decays are helicity suppressed by a factor of mj/m 2 Bc , and hence a precise 
determination of /b c is difficult to make. On the other hand, although the B~ — > tv 
channel is free of the above-mentioned helicity suppression, its observation is possible if we 
have good efficiency. 

When a photon is radiated in addition to the leptons, the helicity suppression is removed 
and a large branching ratio is expected, which is the reason why it makes the investigation 
of the B~ — > iv'j decay much more interesting. 

Note that, the B~ — > £z>7 decay is investigated in the SM within the context of the 
constituent quark model approach ||. But the concept of the "constituent quark mass" is 
itself poorly understood and its relation with QCD is unclear. Therefore the prediction of 
the branching ratio within the context of the constituent quark model approach is strongly 
model dependent. Our aim in this work is to investigate the B~ — > ivj decay in a model 
independent way, namely within the framework of the light cone QCD sum rules (for a 
review on light cone QCD sum rules, see ||). 

The paper is organized as follows. In Section 2 we derive the sum rules for the transition 
form factors, which appears in B~ — > iv^i decay. In Section 3 we present the numerical 
analysis and give a brief discussion about the results. 
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2 Light cone QCD sum rules for transition form fac- 
tors 



The B c — > iv decay at quark level is described by the tree level Feynman diagram, which 
is presented in Fig. 1. The corresponding effective Hamiltonian for the B~ — > iv decay is 



G 

H = -^V cb c'j^ (1 - 75) M'j^ (1 - 75) v 



(1) 



As already noted, the pure leptonic process B~ — > iv [i — e, p) is helicity suppressed. 
If a photon is attached to any of the charged lines, no helicity suppression exists. However if 
a photon is attached to the charged lepton line, it follows from helicity arguments that the 
contribution of such a diagram must be proportional to the lepton mass rri£, and hence it can 
safely be neglected. When a photon is radiated from W^-boson line, the contribution of this 
diagram is also strongly suppressed by a factor m^/ra^, since an extra l/m^y factor comes 
from the second W^-boson propagator. Therefore in B~ — > £P 7 , the main contribution 
should come from the diagrams, where photon is radiated from the initial quark lines. 
Thus the corresponding matrix element for the B~ — > ivj process can be written as 



M = (l(q)\c^ (I-75) b\B c {p + q))£ lf , (1 - 7s ) v 



(2) 



The matrix element (7(9) |c7 M (1 — 75) | B{p + q)) describes the annihilation of the B c mesons 
into the current 07^ (1 — 75) b, with momentum p accompanied by the radiation of a real 
photon with momentum q. This matrix element can be written in terms of two independent, 
gauge invariant structures: 



(7(q)\c^(l-~f 5 )b\B c (p + q)) = e 



m B c 



e* (pq) - (e*p) g M 



(3) 



where e M and q^ are the polarization and the four-momentum of the photon, respectively, 
p is the momentum transfer, g(p 2 ) and f(p 2 ) are the parity conserving and parity violating 
form factors. The main problem is to calculate the form factors g(p 2 ) and f(p 2 ) including 
their momentum dependence, and we will calculate these form factors in the framework of 
light cone QCD sum rules method. 

We start with the following correlator function 



T,(p, q)=i J d 4 xe ipx ( 7 (g) \T {c(x)^ (1 - 75) b(x) 6(0)i 75 c(0)} 







(4) 



We will present the calculation of this correlator in two different ways. In the first 
approach we sandwich T^(p, q) in between the hadronic states with B c meson quantum 
numbers, and hence we get 



T„ip,q) = e 



™ 2 bJb c 
m h + m c 



m B c -{p + qY 



m k 



+ i 



e* (pq) - (e*p) q^ 




(5) 
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where we have used 



(B c \bi l5 c\0) 



m ljB c 

m b + m c 



The second alternative approach, on the other hand, is to calculate the correlator func- 
tion (4) at large Euclidean momentum, where p 2 and (p + q) 2 are both large and negative. 
The Lorentz decomposition of the correlator is 



T»(p, q) = e lxa/3a e* a p p q a T x + % [e* (pq) - (e*p) q^ 



To 



(6) 



The B c meson contains two quarks which interact with the photon only perturbatively. 
This point is essentially different in the B — > lv~f decay in which the photon interacts 
with the quarks both perturbatively and non-perturbatively (see for example (7|). It is our 
aim now to calculate these perturbative contributions. For the invariant structures T\ and 
T2, we write the dispersion relation in variable (p + q) 2 at fixed p 2 . 



7^(1.2) 



ds 



p^ 2 ~>(s,p 2 ) 



(p + qy 



subs, terms . 



(7) 



Here the superscript 1 and 2 corresponds to T x and T 2 respectively. p( 1,2 ) are the spectral 
densities and they are calculated using the method given in || (for applications of this 
method see for example 0, |8|, |9], |i~0f). After lengthy calculations we get for the spectral 



densities 



P 



(i) 



No 
4tt 2 



Qt 



rrib ( In 



1+a-p+X 
l + a-/3- A 



A + m c A 



+Q C 



In 



1 - a + {3 + A 
1 - a + p- A 



- A + m b \ 



(8) 



(2) N c 



4tt 2 (s-p 2 f 



m b Qt 



2m b + p — s J In 



l + a- (3 + \ 
1 + a- f3 - A 



-A (2m 2 - 2m 2 + p 2 {2 - a + (3) - s 



+m c Q b 
+m b Q c 
+m c Q c 



-2m h In 



I + a- f3 + \ 

I I + a-/3- A 



+ A (2m 2 ; - 2m 2 - p 2 {a - /3) + s 



(9) 



2m 2 \n- — Q + / ! + A + A (2m 2 b - 2m 2 - p 2 (a - (3) - s 



1-a+p-X 

1 -a + (3 + A 
a + f3 - A 



(s — p 2 — 2m 2 ^ ln^ — 



A (2m 2 - 2m 2 - p 2 (2 + a - (3) + s) 



where A = \/I + a 2 + (3 2 — 2a — 2/3 — 2af3 and a = m 2 /s, [3 = m 2 /s. Q b and Q c are 
the electric charges of the b and c quarks, respectively and iV c is the color factor and 
5'(s — t) = 45(s — t). Using Eqs. (7), (8) and (9), T\ and T 2 take the following form 
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7i = e 



Nc_ 

47T 2 



s — (p + q) [s — p 2 
l+a-(3 + X 



m b In 



1+a-P-X 



X + m c X 



+ Qc 



(10) 



/ l-a + (3 + X \ 

m c In- — - A + m b X 

\ 1 — a + p — X 1 



T 2 = e 



ds 



4tt 2 

+m c Q b 
+m b Q c 
+m c Q c 



s — (p + q) 2 [s - p 2 ] 2 

1 + a - /3 + X 



(2ml + p 2 — sj In 



-2ml l n 



1 + a - (3 - X 
l + a-(3 + X 
1 + a - (3 - X 



X [2ml ~2m 2 c + p 2 {2 - a + 13) - s 
^^. + X(2ml-2ml-p\a-(3) + s) 



' 2 l-g + /3 + A 

2m c 7^ r 

1 - a + /3 - X 



+ X (2m 2 b ~ 2m l ~ P 2 (® - P) ~ s) 
(s-p 2 - 2m\) ln |l^^-A " A H " ^ ~P 2 ^ + a-P) + s 



a + f3 - X 

To match the results (10) and (11) with the B c meson contribution to the correlator function 
(5), we introduce a new variable u = (m b + m c ) 2 — p 2 / (s — p 2 ) in the dispersion integral 
and invoke the quark-hadron duality in order to subtract continuum contribution, which 
modeled as perturbation contribution starting from some threshold so and finally perform 
the Borel transformation in the variable (p + q) 2 . This operation is necessary to suppress 
higher states and continuum contributions. As a result we arrive at the following sum rules 
for the form factors g(p 2 ) and f{p 2 ): 



9{P 2 ) 



m b + m c N c f 1 du 



ft 



An 2 



e [m% c U-(m b +m c ) 2 + P 2 u]/(M 2 u) 



A U 



(12) 



x 



(Qc - Qb) (m b - m c ) X + Q b m b ln^— ^ — f + ^ + Q c m c In- 1 



1 + a - (3 - X 



I - a + (3 - X 



2 m b + m c N c f 1 



f(p 2 ) = 



du 



x ^m b Q 
+m c Q b 
+m b Q c 
+m c Q c 



f Bc 4vr 2 7a \{m b + m c f _ p z 

1 + a - [3 + X 



3 \m\u- (m b +m c f +p 2 u\ / ( M 2 u) 



(2m 2 h + V 2 



s In 



1 + a - (3 - X 



X (2m 2 b - + p 2 (2 - a + (3) - s) 



l + a - (3 - X 



-2^1n- 1 + a -' ? + A 

1 - 

1 - a + (3 - X 



+ A (2m 2 b - 2m\ - p 2 (a - /?) + s) 



(13) 



1ml In 1 a + P + X + A (2m 2 b ~ ^< -p 2 (a-f3)- s) 



(s-p 2 -2m 2 ) m |_^^ -X(2m 2 -2m 2 c -p 2 (2 + a-(3) + s) j , 
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where in Eqs. (12) and (13) 



_ (m b + m c ) 2 -p 2 u _ _ (m b + m c f - p 2 

s — , u — 1 — u , anct Za — - , 

u Sq — p 2 

and So is the continuum threshold. If we formally set m c — > in Eqs. (12) and (13), the 
resulting expressions are expected to coincide with the form factors corresponding only to 
the perturbative part of the B^ 1 — > £i"y decay. This decay was investigated in 0], and their 
results are identical to ours in the m c — > limit. 

The calculation for the differential decay rate yields 

f = ^V bc \ 2 m\x(l-xf[\g(x)\ 2 + \f(x)\ 2 ] , (14) 

where x = p 2 /m 2 B . The corresponding decay width which follows from the above expression 
is 

„ G 2 q„, „ , r 1 



V bc \ 2 m 3 B [ x(l-xf\\g(x)\ 2 + \f(x)\ 2 \dx. (15) 
Jo 1 J 



96tt 2 ' 

3 Numerical analysis 

In regard to the numerical analysis in the evaluation of the form factors, we have used the 
following set of parameters: m b = 4.7 GeV, m c = 1.4 GeV, mB c = 6.258 GeV |l|, |ll|, [12 
s = 50 GeV 2 and f Bc = 0.35 GeV 0, 0, 0. 



In Fig. 2 we present the dependence of g(p 2 = 0) and f(p 2 = 0) on the Borel parameter 
M 2 . From this figure we see that the best stability is achieved in the range 15 GeV 2 < 
M 2 < 20 GeV 2 for which an uncertainty less than 10% is observed. The analysis for p 2 ^ 
yields similar results as well. 

The sum rules of the type (12) and (13) are expected to work in the region (m?, + m c ) 2 — 
p 2 ~ few GeV 2 , which is smaller than the maximal available p 2 = (nib + rn c ) 2 - To extend 
our results to the whole region of p 2 some extrapolation has been used. The best fits are 
achieved with the following pole formulas: 

v i l—p 2 /mi v 1 l—p 2 /m2 

where 

g(0) = 0.44 ± 0.04 GeV, ml = 43.1 GeV 2 , 

/(0) = 0.21 ± 0.02 GeV, m\ = 48.0 GeV 2 . (17) 

For a more precise determination of the form factors, the perturbative 0(a s ) corrections 
need to be calculated. We will consider this point elsewhere in a future work. 

Fig. 3 depicts the dependence of the differential branching ratio on x. It is observed 
from this figure that the photon spectra (x = 1 — 2E 1 /rriB c ) is practically symmetric. Note 
that the quark model approach also predicts a symmetric spectra for the photon (see Fig. 
2in@). 



5 



Finally, we summarize the numerical results of the branching ratios. Taking \ V c t\ = 0.04 
H, r (B c ) = 0.52 x 10- 12 s pf, we get 



B (B c -> ^7) = 1.0 x 10~ 5 . (18) 

Here we note that the quark model approach yields 

B (B c -> fl/ 7 ) = 5.0 x 10~ 5 . (19) 

If we compare Eq. (18) and (19), it is obvious that the light cone QCD sum rules prediction 
on the branching ratio is approximately five (three) times smaller than the one predicted by 
the quark model approach, if the constituent w-quark mass is taken as 0.35 GeV (0.48 GeV). 

For completeness we present the predictions of the branching ratios for the pure leptonic 
decays f5j: 

B (B c -> nv^) = 6.0 x 10~ 5 , 

B(B c ^eP e ) = 1.4 x 10~ 9 . (20) 

From a comparison of Eqs. (18) and (20), we observe that the branching ratio of the 
radiative leptonic decay is of same order with the corresponding pure leptonic decay (i.e., 
for /i case). 

Since B c — > iuj and B^ —>■ iuj decays have the same final states, it will be interesting 
to know about the relative fraction of the £u , y final states which outcomes from the different 
sources of B c and B ± . 

a) QCD sum rules prediction: 

N B ± 
1 - 0.2 , 



N, 



Be- 



where, for the branching ratio of the B^ — > iwy decay we have used B (B^ — > Iv^j) 
2.0 x 10- 6 0. 

b) The quark model prediction (see M): 



N, 



N 



1.2 (for m u = 0.35 GeV) , 

(21) 

0.7 (for m u = 0.48 GeV) . 



It is quite obvious from the above results that, the two approaches lead to absolutely 
different predictions on the relative fraction of the £i"y final states from B c and B ± . This 
result will be checked in LHC in future experiments. According to the estimations that 
have been made [2-3], the number of B c mesons that are expected to be produced in LHC 
is ~ 2 x 10 8 , and hence the B c — > iv^ decay can easily be detected at LHC which opens the 
way for a real possibility of the experimental investigation of the properties of the B c meson. 
In an investigation of the B c — > iv^i decay, according to our estimation, the background 
contribution from B^ — > is small. 

In conclusion, we have calculated the branching ratio for the B c — > iv^i decay, in the 
framework of the QCD sum rules method within SM and found that B (B c tv^f) ~ 
1.0 x 10~ 5 . 
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Figure Captions 

1. Feynman diagram for the B c — > iv. 

2. Dependence of the form factors g (p 2 = 0) and / (p 2 = 0) on the Borel parameter M 2 . 
In this graph, the continuum threshold s is fixed to the value of 50 GeV 2 . 

3. Dependence of the differential Branching ratio on x for the B c — > iv^i decay. 



7 



Fig. 1 
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